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ABSTRACT 



Context. PSR B0031-07 is well known to exhibit three different modes of drifting sub-pulses (mode A, B and C). It has recently 
been shown that in a multifrequency observation, consisting of 2 700 pulses, all driftmodes were visible at low frequencies, while 
at 4.85 GHz only mode-A drift or non-drifting emission was detected. This suggests that modes A and B are emitted in sub-beams, 
rotating at a fixed distance from the magnetic axis, with the mode-B sub-beams being closer to the magnetic axis than the mode-A 
sub-beams. Diffuse emission between the sub-beams can account for the non-drifting emission. 

Aims. Using the results of an analysis of simultaneous multifrequency observations of PSR B0031-07, we set out to construct a 
geometrical model that includes emission from both sub-beams and diffuse emission and describes the regions of the radio emission 
of PSR B0031-07 at each emission frequency for driftmodes A and B. 

Methods. Based on the vertical spacing between driftbands, we have determined the driftmode of each sequence of drift. To restrict 
the model, we calculated average polarisation and intensity characteristics for each driftmode and at each frequency. 
Results. The model reproduces the observed polarisation and intensity characteristics, suggesting that diffuse emission plays an 
important role in the emission properties of PSR B003 1-07. The model further suggests that the emission heights of this pulsar range 
from a few kilometers to a little over 10 kilometers above the pulsar surface. We also find that the relationships between height and 
frequency of emission that follow from curvature radiation and from plasma-frequency emission could not be used to reproduce the 
observed frequency dependence of the width of the average intensity profiles. 

Key words. Stars: neutron - (Stars:) pulsars: general - (Stars:) pulsars: individual (B0031-07) 



1. Introduction 

Pulsar B0031-07 is well known for its three modes of drift- 
' ing sub-pulses. They are called mode A, B and C and are 
\> characterised by their values fo r of 12, 7 and 4 times 
$_j ■ the pulsar period, respectively (Hu guenin et alj Il970l) . This 
d ' pulsar ha s been thoroughly s t udied at low observ i ng fre- 
quenc i es dHug uenin et al. 1970; Krishnamohanj [l980t I Wright 



198H I Vivekanandl 1 1 995k IVivekanand & Joshi Il997l ll999: 



Joshi & Vivekanandl l2000p. but only rarely at an observing 



frequency above 1 GHz ([Wright & Fowled 1 198 it iKuzmin et alj 



1986; llzvekova et a"illl993l) . Recently. ISmits et alj d2005l) anal- 
ysed simultaneous multifrequency observations from both the 
Westerbork Synthesis Radio Telescope and the Effelsberg Radio 
Telescope and detected all three drift modes at 325 MHz, but 
only detected drift mode A at 4.85 GHz. The pulses that were 
classified as mode B or C at low frequency only showed non- 
drifting emission at high frequency. On the basis of their find- 
ings, they suggest a geometrical model where modes A and B at 
a given frequency are emitted in two concentric rings around 
the magnetic axis with mode B being nested inside mode A. 
This nested configuration is preserved across frequency with 
the higher frequency arising closer to the stellar surface com- 
pared to the lower one, consistent with the well known radius- 
to-frequency mapping operating in pulsars. Due to the rare oc- 
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currence of mode C, they did not attempt to include this drift 
mode in their model. 

Here we analyse new multifrequency observations of 
PSR B0031-07, obtained with the Giant Metrewave Radio 
Telescope, the Westerbork Synthesis Radio Telescope and the 
Effelsberg Radio Telescope simultaneously. In total, the observa- 
tions contain 136 000 pulses spread over 7 different frequencies. 
From these observations we attempt to restrict the geometry of 
this pulsar and create a model that reproduces a great number of 
its observed characteristics. 

In Section|2]we explain how the observations have been ob- 
tained, how the different modes of drift have been determined, 
and which further analyses have been carried out. In Section [3] 
we present our results, which is followed by the modelling of the 
geometry of the emission of PSR B0031-07 in Section |4] The 
discussion and conclusions follow in Sections [5] and [6] respec- 
tively. 

2. Method 

2.1. Definitions 

To describe the observational drift of sub-pulses we use three 
parameters, which are defined as follows: P3 is the spacing, at 
the same pulse phase, between drift bands in units of pulsar pe- 
riods (Pi); this is the "vertical" spacing when the individual ra- 
dio profiles obtained during one stellar rotation are plotted one 
above the other (stacked). P2 is the interval between successive 
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Table 1. List of known param eters of PSR B0031-07. All values 
are from lTavlor etail d!993l) . 



Parameter 


Value 


Pi 


0.94295 s 


P 


4.083 • 10~ 16 


DM 


10.89 pc crrr 3 


>$400 


95mJy 


S 1400 


llmJy 


B sm -f 


6.31 • 10 n G 


E 


1.9 • 10 31 erg/s 



sub-pulses within the same pulse, given in degrees, and A<p, the 
sub-pulse phase drift, is the fraction of pulse period over which 
a sub-pulse drifts in one pulse period, given in °IP\. Note that 
P 2 =PiX A0. 

These parameters are often thought to be associated with 
beams of emission (sub-beams), rotating at a fixed distance 
around the magnetic axis. When each sub-pulse that is seen to 
be drifting in consecutive pulses (which is observed as a drift 
band) is due to emission from an individual sub-beam, then P^ is 
the rotation time of this configuration divided by the number of 
sub-beams. However, when the sub-beams rotate fast enough, it 
becomes possible that multiple sub-beams contribute to one drift 
band. This is called aliasing and is hard to detect. The observed 
P3 is then no longer the rotation time divided by the number of 
sub-beams, but depends on the degree of aliasing. 

2.2. Observations 

Here we present observations of PSR B0031-07, which were 
obtained with the Giant Metrewave Radio Telescope (GMRT), 
the Westerbork Synthesis Radio Telescope (WSRT) and the 
Effelsberg Radio Telescope (EFF). They comprise 7 different 
frequencies. A great part of the observations are simultaneous 
at different frequencies and include 4 hours in which radio emis- 
sion from PSR B0031-07 is observed at 5 frequencies simulta- 
neously. 

The 150, 243, 325, 610 and 1167 MHz observations were 
conducted using the GMRT located in Pune, Indi a. The G MRT 
is a m ulti-element aperture synthesis telescope (Swam p et al.l 
fl99l consisting of 30 antennas distributed over a 2 5 -km diam- 
eter a rea which can be configured as a single dish (Gupt aetal] 
2000). The antennas each have a gain of 0.3K/Jy. The system 
temperature for the different frequencies are T^] = 482 K, T^f 
= 177 K, = 108 K, Tf£ = 92 K and 7^ 67 = 76 K. The sig- 
nals from the antennas at 150, 243, 325 and 610MHz have cir- 
cular polarisation, while at 1 167 MHz the signals are linearly po- 
larised. At any frequency, orthogonally polarised complex volt- 
ages arrive at the sampler from each of the antennas. The volt- 
age signals are subsequently sampled at the Nyquist rate and 
processed through a digital receiver system consisting of a cor- 
relator, the GMRT array combiner (GAC) and a pulsar back-end 
for GAC. The signals selected by the user are added in phase and 
fed to the pulsar back-end. The pulsar back-end computes both 
the auto- and cross-polarised power, which was then recorded at 
a sampling rate of 0.512 ms. We have carried out simultaneous 
observation at 243, 607 and 1 167 MHz by splitting the whole ar- 
ray in to 3 different sub-arrays. We have used a scheme wherein 
the digital sub-array combiner, which does the incoherent addi- 
tion of the multi-channel baseband data from different antennas, 
is programmed to blank the data for a selected set of frequency 
channels for antennas from a given sub-array. This is akin to 



setting non-overlapping filter banks for each sub-array and al- 
lows the 16 MHz bandwidth to be divided between the different 
radio frequency bands of observations, without any overlap be- 
tween signals from the different bands, while still preserving the 
time alignment of the data from the different frequency bands. 
The 243 and 6 10 MHz observations involved polarimetry and 
the observations have been corrected for Faraday rotation, dis- 
persion and parallactic angle. Variations in the rotation measure, 
due to intrinsic fluctuations of the ionosphere during the length 
of one observation causes a deviation in the position angle of 
less than 5° at an observation frequency of 243 MHz. Towards 
higher frequencies, this deviation decreases with A 2 . To correct 
for an instrumental polarisation effect, the polarisation calibrator 
PSR B1929+10 was observed at several parallactic angles, and 
was used to calibrate the data using the technique described by 
Mitr a et all d2005l) . In most cases the corruption due to leakage 
was found to be small, and could be corrected up to an accuracy 
of 5%. 

The WSRT observations were made at a frequency of 
840 MHz with a bandwi dth of 80 MHz and recorded using the 
pulsar back-end, PuMa dVoute et"ail |2002| ). The signals from 
14 telescopes were added with appropriate delays resulting in 
a gain of 1.2K/Jy. The system temperature at this frequency 
is around r, v ,=150K. The Effelsberg observations were made 
at a frequency of 4.85 GHz and a bandwidth of 500 MHz. The 
gain is 1 .5 K/Jy and the system temperature at this frequency is 
T s ys -21 K. The polaris ation calibration is described in chapter 4 
of lHoensbroechtdl999l) . The WSRT observations have been cor- 
rected for Faraday rotation, dispersion and for any instrumental 
p olarisation effects using a pro cedure described in the Appendix 
of lEdwards & Stappersl?2004l) . 

A 50-Hz signal present in the Effelsberg observation has 
been removed by Fourier transforming the entire sequence, re- 
moving the 50 Hz peak and Fourier transforming back. Table [2] 
lists the details of all the observations. All the observations were 
aligned by correlating long sequences of pulses with pulses from 
the observation at 607 MHz that were obtained simultaneously. 
This procedure removes any delay in arrival time between pulses 
at different frequencies due to retardation and aberration. The ac- 
curacy of this alignment is within 1 ms. 

2.3. Finding the drift sequences 

Smit s et ail d2005l) have shown that for at least one of the drift 
modes the value for P3 is the same for both high and low fre- 
quency observations. This is consistent with the concept that the 
emission of pulsars originates from sub-beams of particles that 
rotate around the magnetic axis due to the force-free motion of 
the particles in the strong magnetic field. Here, we have inves- 
tigated the values of Pj for both mode A and mode B over a 
larger frequency r ange, using the same method as described by 
ISmits et all d2005l) . From the observations we chose 1 500 con- 
secutive pulses for which the sub-pulses were visible in at least 
4 frequencies. We then tested for each drift sequence whether 
the value of P3 was the same at each frequency. We can con- 
firm that in these 1 500 pulses, within errors, Pj remains con- 
stant for mode A over the frequencies 243, 607, 840 MHz and 
4.85 GHz. Furthermore, we unexpectedly fo und multiple detec- 
tions of mode-B drift at 4.85 GHz, whereas [Smits et all (12005) 
found no mode-B drift in 2 700 pulses at 4.85 GHz, nor in 5 350 
pulses at an even lower frequency of 1.41 GHz. Still, we could 
not accurately measure P3 at 4.85 GHz when the pulsar was in 
mode B. We did confirm that in this mode P3 remains constant 
over the frequencies 243, 607 and 840 MHz. We then made the 
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Table 2. Details of the observations of PSR B0031-07 



Date 


Start time (UT) 


Telescope 


Freq. 


Time res. 


Bandwidth 


Number of 




(s after midnight) 




(MHz) 


(ms) 


(MHz) 


Pulses 


02-09-2004 


62616 


GMRT 


325 


0.512 


16 


2678 


03-09-2004 


69092 


GMRT 


243 


0.512 


6.25 


14607 


03-09-2004 


69092 


GMRT 


607 


0.512 


9.75 


14607 


03-09-2004 


72150 


WSRT 


840 


0.8192 


80 


16140 


ni no onn/i 
U3-09-2004 


/584/ 


ill 




U.920O 


5UU 


J1U42 


03-09-2004 


80087 


GMRT 


1 167 


1.024 


16 


10145 


04-09-2004 


00554 


GMRT 


243 


1.024 


6.25 


4409 


04-09-2004 


00554 


GMRT 


607 


1.024 


9.75 


4409 


07-09-2004 


69522 


GMRT 


607 


0.512 


9.75 


4718 


07-09-2004 


71441 


WSRT 


840 


0.8192 


80 


16781 


07-09-2004 


73690 


GMRT 


243 


0.512 


6.25 


14394 


24-08-2005 


43212 


GMRT 


157 


0.512 


16 


2229 



assumption that there is indeed one fundamental period associ- 
ated with the vertical spacing between drift bands, such as the 
rotation speed of sub-beams around the magnetic axis, that is 
the same at each frequency at each pulse phase. With this as- 
sumption, we determined the drift mode of each pulse for all 
of the observations. This was achieved by visually finding drift 
sequences at the frequency that showed this drift most clearly. 
Whenever possible, these sequences were then inspected at the 
remaining available frequencies to improve the exact beginning 
and end of the sequences. Finally, Pj was calculated for each 
sequence and classified as either mode A, B, or C. Even though 
the time span of the observations is about 10 hours, we only 
detected 4 short drift sequences with a mode-C drift, none of 
them were detected at 4.85 GHz. Because of the low number of 
pulses which were in a mode-C drift, we limited the analysis to 
sequences of pulses which were in mode-A or mode-B drift. 

2.4. Average polarisation profiles 

Once the drift class of all the sequences was known, all pulses 
in a drift class were averaged together for each frequency sepa- 
rately. This resulted in average mode-A and B polarisation pro- 
files and position angle sweeps at each frequency. We also mea- 
sured the widths of each of these profiles at 10% of the peak 
values. 

2.5. Frequency dependence of the fractional drift intensity 

Smit s et al.l d2005l) have suggested that the mode-A and mode-B 
emission originate from different regions of the magnetosphere. 
This means that the line of sight intersection with the sub-beams 
is different for both drift modes (see their figure 4 and 5). In 
particular, this causes the line of sight to miss the center of 
the ring of mode-B sub-beams at high frequency and only cuts 
through a diffuse component surrounding the sub-beams. As a 
result we do not see any drift bands at high frequency when- 
ever the low-frequency emission reveals a mode-B drift. If the 
line of sight intersection with the sub-beams indeed determines 
how clearly a drift pattern can be observed, then the drift pattern 
can provide information about the geometry of the sub-beams. 
To measure how far the line of sight intersection is away from 
the center of the sub-beams in a drift sequence, we calculated 
the fractional drift intensity as follows. First we took the am- 
plitude at each pulse phase of each pulse in the sequence and 
stacked them on top of each other to produce a 2D time series of 
the same type as shown in Fig. [5] We then calculated the abso- 
lute values of the Fourier transform at each pulsar phase of this 
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Average Power Pulsar phase (Deg) 

Fig.l. Contour plot of the LRFS of PSR B0031-07 at 
0.325 GHz as a function of pulse phase during a sequence of 
mode-B drift. The left panel shows the power spectrum averaged 
over phase. The upper panel shows the power averaged over fre- 
quency. 



2D time series. The resulting Longitude-Resolved Fluctuation 
Spectrum (LRFS) was then averaged over a pulse phase window 
containing the pulsar on pulse, giving a Longitude-Averaged 
Fluctuation Spectrum (LAFS). A contour plot of a LRFS, as well 
as the LAFS, is shown in Fig. [TJ The drift intensity was calcu- 
lated by integrating the LAFS over a small frequency range con- 
taining the reciprocal of the Pj value of the sequence and then 
subtracting the integration of the LAFS over a small frequency 
range containing no periodicities. By dividing the drift intensity 
by the total intensity in the LAFS of the pulsar signal we ob- 
tained the fractional drift intensity. This fractional drift intensity 
was calculated for all sequences of drift at each observed fre- 
quency. When the line of sight starts to miss the center of the 
ring of sub-beams and only cut through diffuse emission, then 
the fractional drift intensity should become low. 

2.6. Width of average drift profiles 

The geometrical model that will be constructed to describe the 
single pulse emission from PSR B0031-07 consists of both 
drifting and non-drifting emission. From the model we can cal- 
culate at each frequency and for each drift mode the width of 



4 



J.M. Smits et al.: The geometry of PSR B0031 



the average profile resulting from the drifting emission only. In 
order to use this width to test how well the model describes the 
observed emission, we needed to determine the width of the av- 
erage profiles resulting from the drifting emission only, in the 
observations. This was done for each frequency as follows. For 
each drift sequence, we determined the average drift profile by 
integrating the power in the LRFS at each pulse phase over a 
small frequency range, containing the reciprocal of the P3 value 
of the sequence. For each drift mode these profiles were aver- 
aged together and their widths were measured. 

2.7. Frequency dependence ofP 2 

The change of P2 with observation frequency reflects the change 
of the size of the emission zone with respect to the total path of 
the line of sight intersection at different heights above the pulsar 
surface. It therefore limits the radius-to-frequency mapping that 
can be used for any model. P2 is often calculated by averaging 
autocorrelations of single pulses, containing more than one sub- 
pulse. For PSR B003 1 -07 this method fails at high frequencies 
where there are no pulses containing more than one sub-pulse. 
We therefore calculated P2, for each drift sequence, as the prod- 
uct of Pj and the sub-pulse phase drift (A0) of that sequence. 
The phase drift was measured by cross-correlating consecutive 
pulses of each drift sequence and averaging the resulting cross- 
correlations. The peaks of the average cross-correlations were 
fitted with a quadratic polynomial, of which the phase shift of 
the maximum was taken as the phase drift. For each drift mode 
and observation frequency the values for P2 were averaged over 
many drift sequences to obtain the frequency dependence of P2 
for each drift mode. The values for P2 are not affected by a pos- 
sible aliasing. 
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Table 3. The 10%- widths of the average profiles from all pulses 
and from the drift modes A and B from the observations of 
PSR B0031-07 at 7 different frequencies. 



Frequency 


width (deg) 


width mode A 


width mode B 


(MHz) 




(deg) 


(deg) 


157 


45.7 ± 1.8 


44 ±5 


46.3 ± 1.4 


243 


42.3 ± 0.6 


39.5 ± 1.5 


42.3 ± 0.6 


325 


40.66 ± 0.15 


39.4 ± 0.3 


40.30 ± 0.12 


607 


39.8 ± 1.1 


32 ±4 


39.3 ± 0.9 


840 


40.9 ± 0.6 


32.2 ± 0.8 


43.0 ± 0.9 


1167 


36 ±3 


25 ±6 


38 ±5 


4850 


33.7 ± 0.8 


23.6 ± 0.6 


32.8 ± 0.9 



Table 4. Average fractional drift intensity of drift mode A and B 
for four different frequencies. 



Frequency 


Fractional 


Fractional 


(MHz) 


drift intensity 


drift intensity 




in mode A 


in mode B 


243 


0.14 ±0.05 


0.20 ±0.05 


607 


0.07 ±0.02 


0.11 ±0.03 


840 


0.07 ±0.02 


0.04 ±0.02 


4850 


0.07 ±0.02 


0.02 ±0.01 



3. Results 

The (polarisation) profiles of PSR B0031-07 for 7 frequencies 
are shown in Fig. [2] No polarisation was recorded for the 157, 
325 and 1167 MHz observations. The 157 and 325 MHz obser- 
vations were not simultaneous with the other observations and 
were aligned to have the center of the profiles at the same pulse 
phase as the center of the 243 and 607 MHz profiles, respec- 
tively. The widths at 10% of the maximum were measured for 
all the intensity profiles. They are shown in TablefJ] The effects 
of dispersion are negligible and have no significant contribution 
to the measured width. The errors quoted here and elsewhere in 
the paper are 1-sigma errors. 

Fig. [3] shows the fractional drift intensity in drift sequences 
for 5 different frequencies over a duration of 4 hours and 15 min- 
utes. Due to interference the drift intensity could not be deter- 
mined in some parts of the observations. These parts and parts 
where no pulses were observed are marked in Fig. [3] as hatched 
areas. The average fractional drift intensities of drift modes A 
and B for four of the frequencies are listed in Table|4] The errors 
were estimated by using the rms of the region in the LAFS that 
does not contain the reciprocal of the P3 value. The 1167 MHz 
observation is not listed because this observation does not con- 
tain enough drift sequences to give a reliable average value. 
The widths of the average drift profiles, i.e. of only the drifting 
component, for each drift mode at each frequency are shown in 
Table|5] At 4.85 GHz there was not enough signal in the mode-B 
profile to measure a width. Table|6]lists the average values of P2 
for drift modes A and B for 7 different observation frequencies. 
Note that the observations at 157 and 325 MHz were not part of 
the simultaneous observations. 



Table 5. The 10%- widths of the average drift profiles for drift 
mode A and B for 7 different frequencies. A '-' indicates that 
there was not enough signal to measure a width. 



Frequency 10%- width of 10% width of 
(MHz) mode A (deg) mode B (deg) 



157 


28 ± 


2 


35.0 ± 


1.5 


243 


27 ± 


5 


35.2 ± 


0.9 


325 


33.4 ± 


0.4 


35.6 ± 


0.4 


607 


27 ± 


3 


30 ± 


4 


840 


12.2 ± 


0.9 


10 ± 


4 


1167 


9 ± 


3 


3.1 ± 


1.1 


4850 


16.4 ± 


0.9 







Table 6. The average values of P2 for drift modes A and B for 7 
different frequencies. A '-' indicates that there was not enough 
signal to measure P2. 



Frequency 


P2 of mode A 


Pi of mode B 


(MHz) 


(deg) 




(deg) 


157 


18.9 ± 


1.2 


19.2 ± 0.7 


243 


24 ± 


4 


18.7 ± 1.6 


325 


19.8 ± 


1.0 


18.7 ± 0.6 


607 


21 ± 


3 


18.7 ± 0.6 


840 


10 ± 


3 


14.0 ± 1.9 


1167 


14 ± 


5 


12 ± 3 


4850 


14 ± 


2 
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Fig. 2. Average polarisation profiles from all the observations. From left to right are the average profiles of all pulses, of pulses that 
are in mode A and of pulses that are in mode B. The solid line is the intensity and the dashed and dotted lines are the linear and 
circular polarisation, respectively. Only the 243, 607, 840, 1 167 and 4850 MHz profiles were obtained simultaneously. All profiles 
are binned to the lowest time resolution of 1.024 ms. The method for alignment is described in Section [2~2l 
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4850 MHz 




9000 
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Pulse number 



1.4xl(F 



1.5xl(T 



1.6x1 



Fig. 3. Fractional drift intensity in drift sequences at 5 different frequencies. The gray columns represent sequences of pulses for 
which a drift intensity could be detected. Light-gray indicates mode-A drift, whereas dark-gray indicates mode-B drift. The hatched 
areas indicate regions during which no pulses were observed, or the quality of the observation was not good enough to determine 
the drift intensity. The lower graph is the continuation of the upper graph. The total duration is 4 hours and 15 minutes. 
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emitting 
field line 



radius of star 



R 

Fig. 4. Sketch of a field line as a function of distance along the 
magnetic axis (z) and distance from the magnetic axis (R) near 
the pulsar, to show the definition of the angle x- Point P on an 
emitting field line is defined by it's polar coordinates r and 0. x 
is defined as the angle between the magnetic axis and the foot 
point of the field line on the surface of the star. 



4. Modelling the geometry 

Here, we set out to find a geometry that can reproduce a great 
number of the observed features of this pulsar, which includes 
the position angle sweep (without the orthogonal mode jumps) 
and the frequency dependences of the width of the average inten- 
sity profile, the width of the average drift profile, the fractional 
drift intensity and P2, for drift modes A and B. 

We first try to limit the geometrical possibilities by fitting the 
widths of the observed average profiles to three different mod- 
els that restrict the relationship between height and frequency 
of emission. All models assume that the intensity of the radio 
emission decreases with distance from the maximally emitting 
ring of field lines as a Gaussian given by 



/ = exp 



2 I X« 



2\ 



(1) 



where / is the intensity, x is the angle between the magnetic 
axis and the foot point on the surface of the star of the field 
line from which emission can be observed (see Fig. |4|, xo is 
the angle between the magnetic axis and the foot point on the 
surface of the star of the maximally emitting field line and Xw is 
the angle between the foot points on the surface of the star of the 
two field lines for which the intensity has dropped to e ~~- times 
the maximum intensity. 

In the first model, the curvature radiation model, we assume 
that the frequency of emission is equal to t he characteristi c fre- 
quency for curvature radiation, given by ( dJacksoni ri999), Eq. 
14.81) 



w = -r — 

2 \r cr 



(2) 



where y is the Lorentz factor of the secondary plasma, c is the 
speed of light and r cr is the radius of cu rvature at the emis - 
sion point. Following the derivation from Gangadhara (2004), 
the dipolar geometry of the magnetic field dictates how the ra- 
dius of curvature is related to the distance between the emission 
point and the center of the star, r, and the polar angle to the mag- 
netic axis, 6, as 



r [5 + 3cos(20)]5 
sin 63 V2[3 +cos(26»)] 



(3) 



For angles smaller than 10° this can be linearized within an error 
of 1.5% to give 



4 

r a = -r/0. 



(4) 



In the second model, the plasma-frequency model, we as- 
sume that the radiation is emitted at the local relativistic plasma 
frequency, given in SI by 



= {yY 



eom e 



(5) 



where n is the sum of electron and positron density and e and m e 
are the election charge and mass. We further assume that there 
is a distinct region around the magnetic axis where shots of par- 
ticles can create a secondary pair plasma with a fixed particle 
density which only depends on altitude. The frequency of oc- 
currence of these shots is assumed to decrease as a Gaussian 
with increasing distance from the maximum intensity field lines. 
Thus, the intensity of the emission falls off with distance from 
the maximally emitting field lines, while the particle density of 
the shots and the frequency of emission remains constant at con- 
stant emission height. Since, along the open field lines, the den- 
sity of relativistically moving secondary pair plasma falls off ap- 
proximately with the cube of the distance to the center of the star, 
in a dipolar field, on account of flux conservation, the frequency 
drops with distance to the power |. This leads to the follow- 
ing relationship between the a ltitude and frequency of emission 
dRuderman & Sutherland!! 19751) : 



(6) 



Final ly, we try a v ersion of an empirical relationship, as 
given by iThorsettl dl99ll) : 



r = h 



wo 



(7) 



which is similar to the plasma-frequency model, but has an ex- 
tta distance parameter ho that allows for a minimum emission 
height. For all three models we simulated the intensity profile at 
7 different frequencies, corresponding to the 7 frequencies of the 
observations, for different values of the parameters of the mod- 
els. These parameters and the limits within which they were var- 
ied, are listed in Table [7] y and a>o were varied logarithmically, 
the other parameters linearly. We then calculated the widths at 
10% of the maximum. From these widths and the widths from 
the mode-A and mode-B intensity profiles from the observations 
we calculated the reduced chi-square for each set of parameters 
for each drift mode. We then minimized these chi-squares to fit 
the values for the parameters. Initial values for a and /3 were 
found by fitting the rotating vector model from Radhakrishnan 
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Table 7. Parameters and their ranges of variation used to fit 
the curvature-radiation model (top), the plasma-frequency model 
(middle) and the empirical model (bottom). 



Parameters for the 


lower limit 


upper limit 


curvature radiation 






model 






a 


0.1° 


6° 


P 


0.1° 


6° 


y 


10 


2000 


Xo 


0.3° 


0.92° 




U.Uo 


U.J / 


Parameters for the 


lower limit 


upper limit 


plasma-frequency 
model 






a 


0.1° 


6° 


P 


0.1° 


6° 


w 


0.01 GHz 


1 000 GHz 


Xo 


0.3° 


0.92° 


Xw 


0.06° 


0.57° 


Parameters for the 


lower limit 


upper limit 


empirical model 






a 


0.1° 


6° 


P 


0.1° 


6° 


ho 


10km 


150 km 




0.01 GHz 


1 000 GHz 


Xo 


0.3° 


0.92° 


Xw 


0.06° 


0.57° 



& Cooke (Radhakrishna n & Cookd[l969l) to the position angle 
from all the observations for which polarisation was recorded. 
As can be seen in the panels of Fig. [2] the position-angle sweep 
is very straight, apart from orthogonal mode jumps which causes 
differences in the position-angle sweep at different frequencies. 
After excluding the pulsar phases that contain an orthogonal 
mode jump, each sweep can be reproduced with the rotating vec- 
tor model from Radhakrishnan & Cooke, when a and are kept 
equal and in the range of 0.1° to 6°. 

Once an initial relationship between the altitude and fre- 
quency of emission was found, we proceeded to include both 
diffuse emission and emission from rotating sub-beams to pro- 
duce single pulses in two different drift modes for a number of 
frequencies. The parameters of this model were then fitted to 
reproduce the position angle sweep and the frequency depen- 
dences of the width of the average intensity profile, the width of 
the average drift profile, the fractional drift intensity and P%, for 
drift m odes A and B. Also, we used the result from lSmits et al.l 
(2005) that shows that for drift mode A at low frequencies, the 
line of sight intersection just penetrates the maximally emitting 
field lines, causing the observed double component in the upper- 
left panel of their Figure 8. 

4. 1 . Results of the Modelling 

For both the curvature-radiation model and the plasma- 
frequency model we did not find any sets of values for the pa- 
rameters to describe the frequency dependence of the profile 
widths correctly. The best set of parameters gave a reduced chi- 
square above 300 and above 15 for the curvature-radiation and 
the plasma-frequency model, respectively. In contrast, we found 
that the empirical model, which includes an extra parameter, de- 
scribed the data with a reduced chi-square ranging from 0.7 to 
3 for many different values for the parameters. Characteristic to 
these sets of parameters were low values for ho and ojq, sug- 



gesting a minimum emission height ranging from 1 to 50 km 
above the stellar surface and a maximum emission height rang- 
ing from 10 to 100 km above the stellar surface, where the low- 
est values for the chi-square always correspond to the lowest 
emission heights. We did not allow emission heights below the 
stellar surface, but such values would give low values for the 
reduced chi-square as well. For comparison, we tried to obtain 
emission heights by modelling the emission for different values 
of the parameters a,xo and;^ w (J3 was kept equal to a) and mak- 
ing no assumptions about the relationship between height and 
frequency of emission. For each frequency, the emission height 
was changed in such a way that the width of the profile from 
the observation would be equal to the width of the profile from 
the model. This always resulted in small emission height differ- 
ences between frequencies, varying from a few kilometers to a 
few tens of kilometers. Also, for each set of parameters it was 
possible to choose values for ho and a>o for which the empirical 
model would fit well to the relationship obtained between height 
and frequency of emission. 

Based on these results we proceeded to make a more com- 
plex model, including two drift modes of sub-pulse emission su- 
perposed on low-intensity diffuse emission, assuming the empir- 
ical model for obtaining the emission height for each frequency. 
The parameters of this model are listed in Table [8] where is 
the radius of the pulsar, D w is the angle between the foot points 
on the surface of the star of the two field lines for which the 
intensity from the diffuse emission has dropped to 60% of the 
maximum intensity, D\ is the intensity of the diffuse emission 
divided by the intensity of the sub-pulse emission at the center 
field line, A^ su b is the number of sub-pulses and P3 is the rota- 
tion time of the carousel divided by iV su b. The superscripts A 
and B refer to the value for that parameter in drift mode A or 
B. Note that we only changed the values of xo and P3 to de- 
scribe a change in drift mode. The rotation time of the carousel 
and the number of sub-beams were derived from the observed 
vertical spacing between drift bands, assuming that the observed 
sub-pulses are not aliased. If we are in fact observing an alias, 
then the rotation time of the carousel as well as the number of 
sub-pulses become less. Table [8] also lists the values of these 
parameters that were obtained after fitting them to describe the 
aforementioned features of the observations. The values for ho 
and a>o in Table|8]lead to emission heights ranging from 2.3 km, 
at 4.85 GHz, to 13.6km at 0.157 GHz above the stellar surface. 
Regardless of the relationship between height and frequency of 
emission, these low emission heights are needed to obtain the 
frequency dependence of P2, which depends on the fractional 
change of emission height. The value for Xn corresponds to the 
last open field line, which for a small value of a is given by 

Fig.|5]shows a gray scale plot of the single pulses from the 
observations on the left and of the single pulses from the model 
on the right, for both the highest and smallest frequencies of the 
simultaneous observations. It demonstrates that the model can 
indeed reproduce both mode A and mode B drifts and also that 
the mode B drift disappears at high frequency. It also shows the 
change in profile width. We further calculated the widths of the 
average intensity profiles and the values for P2 from the model 
for both drift modes at all frequencies. For each frequency, the 
signal-to-noise ratio was adjusted to be the same as the signal-to- 
noise ratio of the data. The widths of both the average intensity 
profiles and the average drift profiles from both the model and 
from the observations for drift mode A and B at each frequency 
are shown in Fig. [6] The average fractional drift intensity from 
both the model and the observations for drift mode A and B at 
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Table 8. Optimised parameter values in a full geometrical model 
for the emission of PSR B0031-07, describing two modes of 
drifting sub-pulses superposed on low-intensity diffuse emis- 
sion, assuming a radius-to-fre quency mapping based on the em- 
pirical model of Ijhorsil {l99l]). 



Values ofP2 



Parameter Value 



Parameter Value 



Pi 

a 
P 

h 

pA 



0.9429 s 
10km 
1.83° 
1.83° 
11km 
1.206 GHz 
13 s 
7s 



N sub 



0.85° 
0.81 ° 
0.13° 
0.22° 
0.001 
9 




Phase (Deg) 



Phase (Deg) 



Phase (Deg) 



Phase (Deg) 



Fig. 5. Gray scale plots of single pulses at two frequencies from 
the simultaneous observations (left panels) and from the model 
(right panels). The upper plots are at 243 MHz and the bottom 
plots are at 4.85 GHz. The first 50 pulses are in drift mode A, 
the following 5 pulses are nulls and the remaining pulses are in 
drift mode B. The pulses in the model were aligned with the 
pulses from the 243-MHz observation. This results in the offset 
between single pulses from the observation and from the model 
at 4.85 GHz. There is some interference visible in the pulses 60 
to 68 of the observation at 4.85 GHz. 

each frequency are shown in Fig. [7] Fig. [8] shows the values of 
P2 from both the model and from the observations at each fre- 
quency. Since the difference in between driftmodes due to 
their different values for^o is far less than the error on the mea- 
sured values for P2, the values of P2 for drift mode A and drift 
mode B were averaged together. Images of the model can be seen 
inFig.El 



24 
22 
20 
18 
16 
14 
12 
10 




0.1 



1 

Frequency (GHz) 



Fig. 8. Value of P2 from the observations (solid) and from the 
model (dashed) at each frequency. Note that the presence of non- 
drifting emission affects the method used to measure ?2- (See 
discussion for more details.) 




Fig. 9. Two close ups of the model of the emission zone of 
PSR B0031-07. The vertical line is the rotation axis of the pul- 
sar. The 7 circles indicate the line of sight trajectories corre- 
sponding to the 7 observed frequencies. The emission zone con- 
sists of 9 sub-beams surrounded by diffuse emission, shown as 
semi-transparent. Both images are to scale. The top image shows 
the location of the sub-beams during mode-A drift. The bottom 
image shows the location of the sub-beams during mode-B drift, 
which lie slightly closer to the magnetic axis. 
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Fig. 6. Widths of the average intensity profiles (upper plots) and average drift intensities (lower plots) from the observations (solid 
line) and from the model (dashed line), for drift modes A (left plots) and B (right plots) at each frequency. 
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Fig. 7. Average fractional drift intensity from the observations (solid line) and from the model (dashed line), for drift modes A (left 
plot) and B (right plot). 
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5. Discussion 

Fig. [2] shows that there is only one polarisation-mode present 
whenever the pulsar is in drift mode A, whereas drift mode 
B shows t wo orthogonal mo des of polarisation, as has been 
reported by ISmits et ail d2005l) . They also report that in their 
325 MHz observations the mode-A profile shows a double com- 
ponent, indicating that the geometry of the emission region dur- 
ing this drift mode is such that over a few degrees of pulse phase, 
one sees emission from a region between the magnetic axis and 
the maximally emitting fiel d lines. This is confirm ed by the aver- 
age intensity profiles from llzvekova et aTl (|T993), where a dou- 
ble component becomes visible at 62 MHz. However, our cur- 
rent observations only show a hint of a double component in 
the mode-A profiles. An explanation for this might be that the 
peak of the mode-A emission zone during the present observa- 
tions lies slightly closer to the magnetic axis than the peak of the 
mode- A emission zone during the observations from Smit s et all 
(2005). This makes it possible that no emission was detected 
from the region between the magnetic axis and the maximally 
emitting field line in the low frequency observations presented 
here. The position angles in Fig.[2]reveal that orthogonal polar- 
isation mode jumps occur at all frequencies. However, in some 
cases the jump is not instantaneous, but rather there is a resolved 
transition, as can be seen at 243 MHz. This is not consistent with 
a mode jump resulting from two competing uncorrelated orthog- 
onal modes of polarisation. Also, the resolved transition, present 
in the 607-MHz observation is in the opposite direction of the 
transition that is visible at other frequencies. This difference is 
not a result from an incorrect polarisation calibration, as it can 
be seen that the underlying position angle sweep is identical to 
that at other frequencies. A brief study of smaller sections of the 
607-MHz observation revealed that the average position angle 
sweep of about 1 000 pulses can contain an orthogonal jump in 
which the position angle either increases within a few degrees of 
pulsar phase, decreases within a few degrees of pulsar phase, or 
changes within a timescale smaller than the timeresolution. 

Further, Fig. [2] shows that the average intensity profiles of 
PSR B0031-07 shift to earlier phase with increa sing frequency. 
This c an also be seen in the average profiles from Izvekov a et al.l 
( 1993). They suggest that the shift can be removed by the adop- 
tion of a higher value of dispersion measure. Since Izvekova et 
al. have determined the dispersion measure by aligning the sub- 
pulses in the center of the profile, which is similar to the method 
of alignment used for the present observations, we briefly check 
if both the shift in their observations and the shift in the cur- 
rent observations could be due to twisted magnetic field lines. 
Since the sub-beams follow the magnetic field lines, the drift 
phase of the sub-pulses will differ between frequencies. Thus 
by "artificially" aligning the sub-pulses between frequencies, a 
shift will be introduced between the average intensity profiles. 
To obtain the observed shift, the magnetic field lines need to be 
twisted around the magnetic axis by approximately 20° between 
the emission altitudes observed at 157 MHz and 4.85 GHz. For 
the present model, this implies a twist of 2° per kilometer alti- 
tude. We can estimate the twist of the field lines due to a current 
of charged particles with a density equal to the Goldreich- Julian 
density, streaming away from the pulsar surface at the speed of 
light and filling the entire polar cap. The pulse phase shift A(f> 
over a height difference Ah is then estimated as 

= _B l _ = 2n_ 

Ah R„ C B Z Pic 



where B^, and B z are the <p and z components of the magnetic 
field, R pc is the radius of the polar cap and P\ is the rotation 
period of the pulsar. For this pulsar this leads to a twist of only 
0.0013° per kilometer altitude. To overcome this difference of a 
factor of 1 000, one needs to assume either unrealistically large 
current densities or much greater emission height differences. 

From Fig.[3]it becomes evident that the mode-B drift is seen 
less frequently towards higher frequencies. Interestingly, even 
though the occurrence decreases with increasing frequency, the 
fractional drift intensity itself of this drift mode does not drop 
until 4.85 GHz and even at this frequency it can still be seen oc- 
casion ally. This is in contrast with the results from Smit s et ail 
(2005) who did not see any mode-B drift sequences in 2700 
pulses at 4.85 GHz, which were observed simultaneously at 
325 MHz, nor in 5 350 pulses at 1.41 GHz, which were not ob- 
served simultaneously at any other frequency. Of course, in the 
absence of a simultaneous observation at low frequency, a mode- 
B drift at high frequency might remain undetected. Also, due 
to the infrequent occurrence of the mode-B drift at 4.85 GHz, 
it is possible that Smits et al. did not detect this drift mode 
at 4.85 GHz because it did not occur during their observation. 
Alternatively, the location of the mode-B emission during the 
present observation might differ from t he location of the m ode- 
B emission during the observations from lSmits et alj (12005). The 
fact that a fractional drift intensity can only be detected sporad- 
ically for mode B at 4.85 GHz, does suggest that the location 
of the mode-B emission can slightly change on timescales of a 
few minutes. Despite these changes over small timescales, the 
average fractional drift intensities that are shown in Table|4]rep- 
resent how the line of sight intersection with the emission region 
changes over frequency. However, it should be noted that resid- 
ual interference and variations in the sub-pulse intensity within a 
drift sequence can affect these values. These effects result in the 
large errors listed in Table |4] Still, the values in Table |4] are con- 
sistent with the suggestion that at all frequencies a clear mode-A 
drift can be detected, whereas the mode-B drift becomes harder 
to detect towards higher frequencies. This can also be seen in 
Table [5] where the widths of the average mode-B drift profiles 
decrease with increasing frequency. The widths of the average 
mode-A drift profiles can also be seen to decrease with increas- 
ing frequency, but even at the highest frequency, the profile is 
still broad enough to allow detection. 

Fig- shows that the model can reproduce the single pulses 
of both drift mode A and B, except for the shift in arrival time 
between the single pulses at low and high frequencies. This shift 
can also be seen in the average profiles of Fig. [2] The change 
of the width of the average intensity profiles with frequency can 
be seen in Fig. [6] It is clear that the widths of the intensity pro- 
files from the model have a steeper spectrum than those widths 
from the observations. This steep spectrum is a consequence of 
the change in height by a factor of 2, which is required to ob- 
tain an observable change in Pi. At 4.85 GHz, the width from 
the model becomes broader again due to the increased inten- 
sity of the broad non-drifting emission with respect to the emis- 
sion associated with the drifting sub-pulses. Fig. [7] reveals that 
the average fractional drift intensity from the model is overall 
much higher than the average fractional drift intensity from the 
observations. This is as expected, since the model does not in- 
clude pulse to pulse intensity fluctuations, which will lower the 
fractional drift intensity. Alternatively, the different spectra of 
the drifting and non-drifting emission might be due to a differ- 
ent emission process rather than a different geometry. However, 
from the results presented here, we feel that it is not necessary to 
assume different emission processes. The model does reproduce 
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the disappearing of drift mode B at 4.85 GHz, as was already 
seen in Fig. [5] Fig. [8] shows the change of P2 over frequency 
from the observations and from the model. At low frequency 
the change of P% from the observation and from the model are 
similar. However at high frequencies, the values of P2 from the 
observations are significantly lower than the values of P2 from 
the model. A steeper spectrum of P2 in the model can be ob- 
tained by making the emission height difference between the 
high and low frequencies larger. However, this would result in 
a steeper spectrum for the width of the average intensity pro- 
files, which does not match the observations. A possible expla- 
nation might be that the method used to measure P2 is influenced 
by the non-drifting emission, which we suggest is present along 
with the drifting sub-pulses. The value of P2 is linearly related to 
the sub-pulse phase drift, which is measured by cross-correlating 
consecutive pulses. In the presence of non-drifting emission, this 
cross-correlation not only peaks at the sub-pulse phase drift, but 
it will also have a broad peak at zero phase shift. When the inten- 
sity of the non-drifting emission becomes strong with respect to 
the intensity of the drifting emission, which happens at high fre- 
quencies, both peaks in the average cross-correlation will merge 
into one peak. This results in a lower value for the sub-pulse 
phase drift and thus also for Pi. At low frequencies we do not 
expect the value for P2 to be influenced by the non-drifting emis- 
sion as it is very weak compared to the drifting emission. And 
indeed, our pr esent values of P2 at low frequencies are similar 
to those fro m Izvekova et alj d 1993b . or those from Fig. 9 from 
Bart eTetaLl £L980), who used auto-correlation functions to di- 
rectly measure P2- Finally, Fig. [9] show the actual geometry and 
emission heights of the emission zone of PSR B0031-07, as 
they result from the model. The lowest and highest circles lie 
at heights of 2.3km and 13.6km from the surface of the star, 
respectively. Fig. [9] shows that the change in the location of the 
sub-beams between mode-A and mode-B drift is very small. 
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6. Conclusions 

We have shown the results from an analysis of a simultane- 
ous multifrequency observation of PSR B0031-07. We have in- 
cluded two non-simultaneous observations to obtain a total of 
7 different frequencies. From these observations we first deter- 
mined the drift mode of each drift sequence at each frequency. 
Contrary to what was expected, we found multiple detections of 
mode-B drift at 4.85 GHz. We then determined the position an- 
gle sweep, the width of the average intensity profile, the width 
of the average drift profile and P2 for each drift mode at each 
frequency. We then tried to fit three emission models to the 
widths of the average intensity profiles to ultimately find a model 
that describes all the observed features. We found that the mod- 
els based on curvature radiation and plasma-frequency emission 
could not reproduce the frequency dependence of the widths. 
However, the model based on an empirical relationship between 
the height and frequency of emission, which includes an ex- 
tra parameter, could reproduce the frequency dependence of the 
widths for many different values of the parameters. We therefore 
improved the latter model to reproduce two drift modes of the 
single pulse emission. 

We can summarize the features of the geometrical model of 
PSR B0031-07 that is presented here as follows. 

- The model reproduces to great extent the position angle 
sweep (without the orthogonal mode jumps) and the fre- 
quency dependences of the width of the average intensity 
profiles, the width of the average drift profiles, the fractional 
drift intensity and P2, for drift modes A and B of the single 
pulses of PSR B0031-07. The largest deviations are found 
in the average intensity profiles. 

- The emission heights are very low. The high frequency emis- 
sion comes from a region just above the surface of the star. 
The low frequency emission comes from a region about 10 
kilometers higher than the high frequency emission. 

- The parameters a and /3 are approximately the same and de- 
pending on the actual emission height, around 2° to 4°. 

- The emission is centered around, or close to the last open 
field lines. 

- The emission from drift mode B comes from a region just 
slightly closer to the magnetic axis than the emission from 
drift mode A. 

- Along with the drifting sub-pulses there is non-drifting emis- 
sion in the single pulses that becomes more significant to- 
wards higher frequencies. This non-drifting emission might 
provide a hint as to why not all pulsars show clear drifting 
sub-pulses. 

- Assuming that the observed drift speeds of the sub-pulses are 
not aliased, the number of sub-beams is around 9. 

The model results in very low emission altitudes, ranging 
from 2.3 to 13.6km above the surface of the star. This is in 
strong contrast with other emission heights that have been 
meas ured for pulsars, which are typically some 10 to 1000 km 
(e.g. von Hoensbroech & Xilouris 1997; Gangadhara & Gupta 



2001; Mi tra & Rankin! [2002; MitrajfclJ |2004|). However, by 
improving the model, the emission heights can become higher. 
The model might be improved by assuming a more realistic 
particle-density distribution near the polar cap, resulting in vari- 
able emission heights at a fixed emission frequency. This can 
also improve the plasma-frequency model. Also, the emission 
altitude of mode A and mode B need not be the same, as has 
been assumed here. 
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